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1 Introduction 



It has been recently observed [|l], ^ that most of nonhnear integrable equations, not only 
two-dimensional ones, possess a special discrete symmetry, that appears very useful in many 
respects. In Ref. we demonstrated that this symmetry considered as a transformation 
of the phase space of the system is a canonical transformation. This fact was one of the 
reasons to call the considered transformations the canonical symmetry. In the next paper Q 
we considered the behaviour of the densities of local conservation laws under the canonical 
symmetry for the case of the nonlinear Schrodinger equation. These densities appeared to 
change under the action of the canonical symmetry by the total space derivatives. We have 
noted in that this observation does not give us a constructive method to find the local 
conservation laws for the system under consideration. 

In the present paper we continue the investigation of the properties of the canonical 
symmetry from the Hamiltonian point of view. In particular, we show how the canonical 
symmetry can be used to look for different Hamiltonian operators relevant to the system. It is 
known that if two Hamiltonian operators, that can be used to write the considered evolution 
equations as the Hamilton equations, form a Hamiltonian pair, then we can find the whole 
hierarchy of the involutive conservation laws, starting from the given one 0. It appears 
that the use of the canonical symmetry leads us to a set of the Hamiltonian operators, any 
two of which form a Hamiltonian pair. Thus, the canonical symmetry not only gives a new 
characterization of the local conservation laws but also helps us to find them. 

In this paper we use the notations and conventions of our paper Thus, we consider 
the case of two independent variables t, x and A dependent variables m", a = 1, . . . , A. We 
denote by A the space of the functions of the variables x, m", and the derivatives of m° over 
X up to some finite order. For a general function of such type we use the following notation 

/M = /(a;,M,M, m'), (1.1) 

(fc) (fc) 

where for any k u is the set formed by the k-th derivatives m " of the dependent variables u"" 
over X. We assume also that u'" = u"'. Recall that we call the Hamilton equations evolution 
equations, that can be written in the form 

«-J-M^. (1.2) 

where J is a Hamiltonian operator In the present paper we consider Hamiltonian 

operators that include inverse powers of the operator of the total derivative over x. To 
deal with such operators we have to use a generalization of the usual binomial coefficients. 
This generalization is reminded of in section 2. In section 3 we give the necessary facts on 
integro-differential operators. Section 4 is devoted to the consideration of the invariance 
of a matrix integro-differential operator under the action of a differential transformation. 
This condition for the case of a Hamiltonian operator is equivalent to the requirement of 
the canonicity of the transformation. In section 5 we consider a partial case of the second 
order transformations. Using as examples the nonlinear Schrodinger equation and modified 
nonlinear Schrodinger equation we find all matrix skew-symmetric operators which satisfy 
the invariance condition for the corresponding canonical symmetry. It appears that they 
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are actually the Hamiltonian operators, forming the well known hierarchies of Hamiltonian 
operators for the nonlinear equations under consideration. 



2 Binomial Coefficients 



Here we define the binomial coefficients y^j for arbitrary integers n and k, such that n > k. 
Recall that for n,k > and n > k we have 

[k J k\{n-k)V ^ ■ ' 



Let us write the well known relation 



(i+^r = E(^W'- (2-2) 



k=0 



It is convenient to put for < 

\k, 



0. (2.3) 



Hence, we can rewrite Eq. (2^) in the form 



{l + xT= E {\x\ (2.4) 

k=—oo \ / 

We use this formula as the definition of the binomial coefficient for n < 0. In particular, for 
n = — 1 we have 

(l + x)-i= E (-l)-'=-V, (2.5) 

k=—oo 

hence, for k < —1 

^"^^ =(-1)-^-^ (2.6) 



k 



For any n < we can write 
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Differentiating Eq. (p.5|), we get 



\k J \ —n — 1 

From this relation we have subsequently 



( n \ ( n \ nin — 1) , , 



It is clear that for any A; > 

( n \ n{n — 1) ... (n — A; + 1) 
\n — kj k\ 

Note, that equality ( |2.1UD is valid also for n > 0. 



(2.10) 
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3 Integro— Differential Operators 

In this section we shall give necessary facts on integro-differential operators. A more detailed 
exposition can be found, for example, in . An integro-differential operator is a formal series 
of the form 

N(a) 

a= J2 (3.1) 

k=—oo 

where Ofc G A. We suppose that aN{a) 7^ 0, and say that the number N{a) is the order of the 
operator a. 

The multiplication of integro-differential operators is defined as follows. The product of 
two integro-differential operators a and h is the integro-differential operator 

N{c) 

c = ab= CkD\ (3.2) 

fc=— oo 



where N{c) = N{a) + N{b), and 

Nia) 



I / J\ 

l=k-N(b) m=k-N(b) V^J 



k—m ) 



E E [\^^]ak-n.D^-{h). (3.3) 



l=k-N{a) m=k-N{a) 



This relation can be formally obtained by the multiplication of the corresponding formal 
series and by the appropriate resummation. The direct calculation shows that thus defined 
product is associative. 

The transposed operator of the integro-differential operator a is the operator 

N(a) 

«^ = E i^^)kD', (3.4) 

fc=— oo 

where 

N{a) /,\ 

{a^)k^j:i-iy^f^jD^-Hai). (3.5) 

In particular, we have 

D'^ = -D. (3.6) 
We shall need below the following relation 

N{b) l+N(a) / \ 

iabnk= E E i-iyr]a^-iD-Hh) 

l=k-N(a) m=k \ / 



N{a) 
N{a) l+N{b) 

E i-^r'TuhiD-'-'ib^-i)- (3-7) 



m 

(-1)'— [ 

l=k-N(b) m=k 
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An integro-differential operator a is said to be invertible if there exists an integro- 
differential operator b such that 

ba = 1. (3.8; 

If it is the case we write b = a~^. It is clear that N{a~^) = —N{a). 

Suppose that an operator a is invertible, then from Eq. (|3.3|) for A; = we get 



(a ^)~N{a)0'N{a) = 1- (3.9) 

Hence, the operator a is invertible only if 

«7V(a) 7^ 0, (3.10) 

and in this case 



(a-i)_;v(a) = (3.11) 



From Eq. ( p.3|) for < we have 



-N{a) 



{a-')k-Nia)aN{a)+ E E ( ^ ) (a"')i/^'"™(afc-™) = 0. (3.12) 

l=k-N{a)+l m=k-N{a) ' 



We can consider Eq. ( p.l2[ ) as a recursive relation for the determination of the quantities 



(a~ )fc, and conclude from it that the operator a is invertible if and only if Eq. ( |3.10| ) is 
satisfied. 

The matrix integro-differential operator A is a matrix with the matrix elements being 
integro-differential operators. The transposition operation can be naturally extended to the 
case of such operators. It is clear also that a matrix integro-differential operator 

N{A) 

A= E ^fe^' (3-13) 

k=—oo 

is invertible if and only if the matrix An{a) is invertible. 

4 Invariance Condition 

Recall [^] that the invertible differential transformation 

= = ^-(x, u, u, %) (4.1) 



is a symmetry of Hamilton equations (|1.2|) if 

ip'[u]J[u]ip'^[u] = J[^[u]] (4.2) 

and 

h[ip[u]] - h[u] e Ker 6/6u. (4.3) 
In Eq. ( |4.2| ) if' denotes the Frechet derivative of the differential transformation ( ^.1| ) 
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In papers ^ we have shown that the canonical symmetry of a nonhnear system is a 
symmetry of the corresponding Hamilton equations. It is known that there are different ways 
to write integrable equations in the Hamiltonian form. In this the role of the corresponding 
Hamiltonians is played by different conserved quantities. We proved in 0], for the case of 
nonlinear Schrodinger equation, that the densities of local conservation laws change under the 
action of the canonical symmetry by total space derivatives. In other words, if we consider 
the density of a conservation law as the density of the corresponding Hamiltonian, then we 
see that it satisfies Eq. ( |4.3| ). It is natural to suppose that the corresponding Hamiltonian 
operator satisfies Eq. ( [4.2[ ), which we call the invariance condition. We shall not demonstrate 
that directly, but consider a more general problem: what the Hamiltonian operators are that 



satisfy Eq. (|4.2| ) for given differential transformation (4.1). 

In this paper we consider only the case of two dependent variables and use the notation 
q = M^, r = u^. In this case the operator J can be written as 




(4.4) 



The condition J = —J^ gives 



thus, we can write 



= -a, (f = -d, (4.5) 
J = -b, = -c, (4.6) 



(4.7) 




The operator (f' in the case under consideration may be represented as 

and invariance condition ( |4.2| ) is equivalent to the relations 

a =aaa^ - (3b^a^ + ab(3^ + (3d(3^, (4.9) 

b =aa-f'^ - (3b^-f^ + ab6'^ + (3d6'^, (4.10) 

- F= 7aa^ - + 76/?^ + 6d(3^, (4.11) 

d =^a-f'^ -6b^-f^ + -fb6^ + 6d6^, (4.12) 



where a[u] = a[ip[u]], etc. Note that Eq. (|4.11|) can be obtained from Eq. ( |4.1CI| ) via the 



transposition. Thus, we have only three independent relations. In all the cases that we 
consider in this paper a = and 7 are simply functions. Therefore, we can rewrite system 
(^-(TO) in the form 



a- (3d(3 = 0, (4.13) 
5d + + 76 = 0, (4.14) 

5/3-^6 + 765^ + 707-^ = 0. (4.15) 
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From this system we conclude, in particular, that 

N{b) = N{a) + N{6), N{d) = N{a). 

Representing a and d as 

N{a) N{d) 

k=—oo k=—oo 



(4.16) 



(4.17) 



and using Eq. ( |3.3|) , we get from Eq. ([4.13|) the following equalities 



m=k \ / 



< N{a). 



(4.18) 



In the same way, using Eqs. ( p.3|) and (|3.5| ), we write Eq. ( [4.14| ) in the component form 
{k < N{h)): 



N(5) I 

E E 

l=k-N{d) m=k-N{d) 



I 



(5;D'-"^(4-„^) 



N{b) 

+ E(-i: 

m=k 



m ' 
, k I 



(4.19) 



Using now Eq. (3^), we get from Eq. ( 4.15 ) for N{a) < k < N{b) + N{S) the relations 



A^(<5) I 

E E 

l=k-N{b)m=k-N{b) 



\S,D'-'"{p-X-n,) 



+ (-i)'('"+,'^"'>)7n-"'"-'(*)*™.«,.,-,)=o. 



(4.20) 



while for k < N[a) we have 

N{5) I 

E E I 

l=k-N{b) m=k-N{b) 



'm+N{b)-l 



N(a) 



E (l]lD^'\^)a^-dk = Q. 

m=k 



(4.21) 



5 Invariance Condition for Second Order Transforma- 
tions 
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5.1 General Consideration 



In this section we consider invariance condition ( [4.13| )- (|4.15| ) for the case N{6) = 2, i. e. for 
the case when 

S = 52D^ + 5iD + 5o. (5.1) 



(0) 



First, suppose that there exists a Hamiltonian operator J of form (|4.7| ) with the ma- 
trix elements being functions. In other words, we suppose that there exists a Hamiltonian 
operator of the order 0, having a special form. From Eq. ( [4.16|) it follows that in this case 
a = d = 0. Denote the function, playing the role of the operator b, by g, and consider 
Eqs. ( [4.13|) -( ^.15|) . Eq. ([4.13|) is satisfied by definition. In the case under consideration 
Eq. ( [4.14| ) takes the form 

9 + (3i9 = 0. (5.2) 

Note, that in all known cases the canonical symmetry has the following property. The 
function f[u], satisfying the relation 



fMu]] = f[u] 



(5.3) 



is simply a constant. From this property it follows that if there exists a function g, satisfying 
Eq. ( [5.2| ), then this function is unique up to a constant factor. 
Eq. (|4.15| ) in our case is equivalent to 

h9 = 195'^- (5.4) 

It is easy to show that the function g satisfies Eq. ( |5.4|) if and only if it satisfies the relation 

D{82)ig - hD{^g) - d^g = 0, (5.5) 

that can be written as 

(5.6) 



\19, 



19' 



Proceed now to the investigation of invariance condition ( [4.13D -( ^.15[) , for the case when 
the order of the operator b is an arbitrary integer N. Denote the corresponding operator by 



(N) 



J . It is a priory unknown, that this operator exists. On the other hand, if it exists then it 



(fc) 



is defined at least up to an arbitrary linear combination of the operators J with k < N. 



(JV) 



Let us suppose that the operator J exists. From Eq. (|4.20| ) in the case k = N + 2 we 
get 

&7V + /?7&7V = 0. (5.7) 
Thus for any the function fo^r is determined up to a constant factor and is proportional to 



(5. 



the function g. Thus, taking into account Eq. (|5.5| ) we see that 

D{52)ibN - 62D{-fbN) - hlbN = 0. 
For k = N + 1 from Eq. ( ^.20| ) we get the following relation 

^2(^^-1 - (3lbN-i) + {{N + 2)D{62) - 5^)l3ibM^ 

+ 5ibM + 252(3 D{(3-^bN) = 0. 



(5.9) 
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Using Eqs. (|5.7| ) and (|5.8|), we can rewrite this equality in the form 

52(pN-i + hhN^i) + N(5-iD{52)hM = 0. (5.10) 
It is convenient to introduce the functions 

fk = bk/bN. (5.11) 
For the function fN-2 from Eq. (|5.10|) we have 



fN-i-fN^i = N^. (5.12) 

Hence the function /n-i, if it exists, is defined up to an additive constant. The function 
bN-i, in turn, is defined up to a constant multiphed by the function bN- Recall that the 
function fo^v is proportional to the function g. It is clear, that the operator J{N) is defined up 
to a linear combination of the operators J{k) with k < N. It is the reason of the ambiguity 
that we have come across. 

After some rather long calculations from Eq. ( [4.20| ) for the case k = N we get the 
following equality 



S2{bN-2 + hbN-2) - 262f3^D bN - N52hD 



N J \^2^ 
1 



+ (AT - l)f5^D{52)bN-i + ^N{N - l)f3-fD\62)bN = 0, (5.13) 
that can be written in terms of the functions fk as 

fN-2 - fN-2 = - 2D{fN-i) - ND j 

, _ , l^^^ (5,4) 

t)2 2 02 

We see that the function bN-2 is defined ambiguously, and the reason of this ambiguity is 
the same as for the function &iv_i (see above). 
Consider now Eq. ( [4.19[ ). For k = N we get 

rfiv-2 = ((-1)^-1)^. (5.15) 



Thus for an odd we have 



and for an even A^ 



dN-2 = 7 — , (5.16) 



dN-2 = 0. (5.17) 



Note, that we have not actually used the fact that the operator d must be skew-symmetric, 
but Eq. (|5.17| ) is already in agreement with it. 
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We get further 



d^_3 = -{N-2)d(^]= ^i—^D{dr,-2) (5.18) 



for an odd A^, and 



dN-, = - - -^^^ (5.19) 



S2 SI 

for an even N. Eq. ( ^.181) is again in agreement with the skew-symmetricity of the operator 
d. 

At last we get 



dN-, = - + ' ':r + (iv + i)D . ^ 



2 

+ - i(iV^ - 3iV + A)D' ( ^ ) (5.20) 



for an odd A^, and 



A^ - 3 

dN-4 = -^—DidN-3) (5.21) 



for an even A^. 

The functions can be determined from Eq. ( 4.18|) . 

Let us make here some conclusions from our consideration of the invariance condition. 

(AT) 

The ambiguity in the definition of the operator J can be described as follows. Let for any k 

the operator J be some solution of the invariance condition. A general operator, satisfying 
the invariance condition, and having the order A^, is given by the formula 

J= ^kJ, (5.22) 
fc=— 00 

where Uk are some constants. 

(0) (1) 

If we have two solutions, for example J and J, of the invariance condition, and the 

(0) 

operator J is invertible, then the operator 

(2) (1)(0) (1) 

J=JJ-^J (5.23) 

is of the order two and also satisfies the invariance condition. Analogously, we get the 

(1) 

operators, having any positive order. If the operator J is invertible, then we can construct 
an operator of any negative order. For example, the operator of the order —1 has the form 

(-1) (0){1) ,{0) 

J = JJ-^J. (5.24) 
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5.2 Nonlinear Schrodinger Equation 

As a first concrete example, we consider the nonlinear Schrodinger equation. In fact to define 
the canonical symmetry we should consider the following complex extension of it [Q: 

iq + q" - 2erq^ = 0, ir - r" + 2egr^ = 0, (5.25) 

where q and r are arbitrary complex functions of the variables x and t, e is the coupling 
constant. In Eq. ( p.25| ) and below dot and prime mean the partial derivative over t and x, 
respectively. The canonical symmetry for this system has the form ^, ^ 

q = — , r = er^q — r" H . (5.26) 

er r 

Thus, we have 

13 = -\ 7 = er^, (5.27) 

S = -D^ + 2^D + 2erq-^. (5.28) 

Equations (|5.25| ) can be written in Hamiltonian form (|1.2|) if we choose 

(0) / _^ \ 

J=J=(° (6.29) 

and 

h[q, r] = r'q' + er^g^. (5.30) 

(0) 

It is not difficult to show that the operator J satisfies the invariance condition p, Hence, 
in this case the function g is a constant function. 

Recall 10, ^ that for the conserved quantity rq we have 

rq = rq — D i — ) . (5.31) 
\erj 

Consider now the case = 1. From Eq. (|5.7|) we get 

bi-bi = 0, (5.32) 

(1) 

hence, the function bi is a constant function. It is convenient to normalize the operator J 
by the condition 

bi = -1. (5.33) 

From Eq. ( ^.10|) we have 

bo-bo = 0. (5.34) 

Let us put 

bo = 0. (5.35) 
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Eq. ( ^.13 ) is more nontrivial: 



b-i - + 2D 



0. 



(5.36) 



Comparing this relation with Eq. ( |5.31| ), we see that we can choose 

6_i = 2erq. 
From Eqs. ( ^:T6D , ( ^ISj) and ([CTp we find 

= -2er^ d_2 = 2err', d.3 = -2err", 

while Eq. ( |4.18| ) gives 

a_i = -2eg^, a_2 = 2eqq', a_3 = -2eqq" . 



(5.37) 



(5.38) 



(5.39) 



Taking into account the expressions obtained for 6^, dk and we can suppose that the 



(1) 



operator J has the form 



(1) 
J 



-2eqD-^q -D + 2eqD-^r 
-D + 2erD-^q -2erD-^r 



(5.40) 



A direct calculation shows that the operator, given by Eq. ( p.40| ) satisfies the invariance 
condition. Thus, we get the well known second Hamiltonian operator for the nonlinear 
Schrodinger equation H, 0] as a solution of the invariance condition. 



5.3 Modified Nonlinear Schrodinger Equation 

We consider here the following complex extension 0, 0, 0] of the modified nonlinear 



Schrodinger equation |10 



iq + q" — 2ie{rq)q' = 0, ir — r" — 2ie{rq)r' = 0. 
The canonical symmetry for this system is the transformation of the form [|l], ^ Q 

^ 1 



er 



r = erq + i\ r — 



Thus, in this case 



P = - 



7 = er 



i-D^ + i[l + — ]D + 2erq - i—. 



(5.41) 
(5.42) 

(5.43) 
(5.44) 



Equations (|5.41|) can be written in Hamiltonian form (|1 . 2|) for the Hamiltonian operator 
J, given by Eq. ( |5.29| ) and 



r] = r'q + "^('^ QQ' ~ Q ^f')- 



(5.45) 
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(0) 

It is not difficult to show that the operator J satisfies again the corresponding invariance 
condition [Q], and the function g is again a constant function. 

Write now the formulae, describing the behaviour of the densities of the lowest conserva- 
tion laws under the action of the canonical symmetry: 

(5.46) 

irq = —irq' + D ( irq H — - ] . (5-47) 



~q = rq-\ — D I In — 



er 



(1) 



Proceed now to the construction of the operator J . The relation, determining the func- 
tion 6i, is again (|5.32| ), and we can put 



h, = -1. 



From Eq. ( ^.12|) we find 



/o-/o-/^ In- 



0. 



Taking into account Eq. ( ^.46|) , we see that we can put 

ho = ierq. 

Eq. ( ^liD takes in our case the form 



/_i - /_i - eD 2irq + - + — 
V er er 



0. 



Using Eq. ( |5.47|) , we find the following expression for 

6_i = ie{rq' — r'q). 

As in the case of the nonlinear Schrodinger equation we obtain 

d_i = 2ierr', d-2 = —i^ir'"^ + rr"), = ie{rr"' + r'r"); 

o-i = -2ieqq', a_2 = ie(g'^ + qq"), a_3 = -ie{qq"' + q'q") 

(1) 



Thus, we can try as the operator J the operator of form ( ^.7| ) with 

a = -ie{qD-^q' + q'D-^q), 

b = —D + ierq + ie{rD^^q' — r'D~^q), 

d = ieirD'^r' + r'D'^r). 



(5.48) 
(5.49) 

(5.50) 
(5.51) 
(5.52) 



(5.53) 
(5.54) 



(5.55) 
(5.56) 
(5.57) 



A direct check shows that this operator satisfies the invariance condition. 
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6 Conclusion 



In this paper we have shown how the canonical symmetry can be used for the construction of 
Hamiltonian operators relevant to the system under consideration. Note, that we have used 
only the invariance condition, and discovered that the matrix skew-symmetric operators, 
satisfying this condition, are the required Hamiltonian operators. These operators define the 
Poisson brackets, satisfying Jacobi identity, and any two of them form a Hamiltonian pair. 
It is clear that for an arbitrary differential transformation this will not be the case. An open 
problem here is to formulate conditions on the transformation that provide such results. 
Apparently this problem is directly connected to the problem we stated in 0], that was: 
what are the conditions providing the involutivity of the quasi-invariants of a differential 
transformation? 
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